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We show that for an m-partite quantum system, there is a ball of radius 2~' m/ ' 2 ~ 1 - t in Frobenius 
norm, centered at the identity matrix, of separable (unentangled) positive semidefinite matrices. 
This can be used to derive an e below which mixtures of e of any density matrix with 1 — e of the 
maximally mixed state will be separable. The e thus obtained is exponentially better (in the number 
of systems) than existing results. This gives a number of qubits below which NMR with standard 
pseudopure-state preparation techniques can access only unentangled states; with parameters real- 
istic for current experiments, this is 23 qubits (compared to 13 qubits via earlier results). A ball of 
radius 1 is obtained for multipartite states separable over the reals. 
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I. INTRODUCTION 

Entanglement is an important quantum resource, use- 
ful in quantum computation, cryptography, and commu- 
nication protocols. Entangled quantum states are those 
that cannot be expressed as a mixture of product states. 
That is, if p is an entangled state of m systems (an "Tri- 
partite state"), there is no way to choose probabilities Pi 
and states p\ , p™ for systems 1 through m, such that 

t 

In this paper, we provide a simple geometric condition 
sufficient to guarantee separability (non-entanglement) 
of an m-partite state: that the state is proportional 
to the identity matrix plus a Hcrmitian perturbation 
A whose Frobenius norm (2-norm) is no greater than 
2-(m/2-i) This is exponentially better than the best 
previous bounds we are aware of mm it can be used to 
obtain balls of normalized separable states. Because the 
set of separable (unentangled) density matrices is con- 
vex, and the size of the largest ball that fits inside it 
(as well as the smallest ball that covers it) is important 
to complexity-theoretic questions involving convex sets, 
we expect the result to have applications in complexity 
questions about entanglement, such as the complexity of 
deciding whether or not a multipartite state is entangled. 
Equally importantly, it can help determine whether or 
not entanglement is present in interesting theoretical and 
experimental situations. For example, though the util- 
ity of the criterion is emphatically not restricted to such 
states, it gives us a bound on the "polarization" e below 
which "pseudopure" states of the form (1 — e)I/d + en, 
with 7T pure, I the identity operator on a multipartite 
state-space of overall dimension d are separable. Applied 
to a standard nuclear magnetic resonance (NMR) quan- 
tum information processing (QIP) protocol, pseudopure- 
state liquid NMR, at 300 Kelvin and all Tesla external 
field, it tells us that unless we have 23 or more nuclear- 
spin qubits, only separable states can beproduced (com- 
pared to 13 qubits using the bound in |l|). 



The methods of [l|,|2j are very different from ours: they 
expand the density matrix in an overcomplete basis of 
pure states, and find conditions guaranteeing positivity 
of all coefficients of the expansion, thus giving an explicit 
decomposition of the form (JJ. In contrast, the methods 
we use have a nonconstructive flavor: although they es- 
tablish that any m-partite (unnormalized) density matrix 
within a distance l/2( m / 2_1 ) is separable, they do not 
provide an explicit "separable representation" (JTJ of it. 
Our methods use general concepts of matrix theory and 
convex analysis in terms of which the problem is naturally 
formulated, and involve only short and elementary calcu- 
lations. This paper is a natural sequel to Q, in the sense 
that almost the same mathematics is used. The main 
(though quite simple) novelty here is a generalization of 
separability, so-called (Ci(g>C2<B>-"<B>C m )-separability , 
where Ci are matrix cones (see Definition 2 below). This 
generalization arises naturally in extending the bipartite 
result of to m-partite systems with m > 3. A reader 
comfortable with the technique used in Q should have 
no extra problems in understanding this paper. 

II. NOTATION AND MATHEMATICAL 
PRELIMINARIES 

The mathematical notion of a "regular" positive cone 
(which we will just call cone) is basic in quantum infor- 
mation science, especially in the study of entanglement. 
This is so because the unnormalized quantum states, the 
unnormalized separable states of a multipartite quantum 
system, the completely positive maps, the positive maps, 
and many other sets of interest form such cones. (Appro- 
priate normalization-like conditions, such as unit trace 
for states, or trace-preservation or trace-nonincrease for 
maps, are usually just additional linear equalities or in- 
equalities.) In this section, we review regular positive 
cones and related notions; background and preliminar- 
ies specific to the separable cones (of unentangled states) 
appear at the beginning of the next section. 

Definition 1 A positive cone is a subset K of a real 
vector space V closed under multiplication by positive 
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scalars. It is called regular if it is (a) convex (equiva- 
lently, closed under addition: K + K = K ), (b) generat- 
ing (K — K = V , equivalently K linearly generates V ,) 
(c) pointed (K n —K = 0, so that it contains no non- 
null subspace ofV), and (d) topologically closed (in the 
Euclidean metric topology, for finite dimension). 

Such a positive cone induces a partial order >^ on 
V, denned by x >k y ■— x — y G K. It is "linear- 
compatible": inequalities can be added, and multiplied 
by positive scalars. A set S is said to generate a cone K 
if K is the set of positive linear combinations of elements 
of S. The topological closure condition guarantees that 
such a cone is generated (via addition) by its extreme 
rays. These are sets R x := {Xx : X > 0} such that 
no y G R x can be written as a convex combination of 
elements of C that arc not in R x . We will not make 
much use of closure and extremality, but at some points 
we use the fact that positive semidefinite (PSD) matrices 
can be written as convex combinations of rank-one PSD 
matrices (these being the members of the extreme rays 
of the cone of PSD matrices). 

Duality is often a useful tool when dealing with cones. 
The dual vector space V* for real V is the space of linear 
functions ( "functionals" ) from V to R; the dual cone 
C* C V* of the cone C C V is the set of such lin- 
ear functionals which are nonnegative on C. For finite 
dimensional vector spaces, V* is isomorphic to V as a 
vector space (they have the same dimension). However, 
this isomorphism is not canonical; any nonsingular linear 
maps from V onto V* does the job. For a slight improve- 
ment in clarity below, for a G V* , x G V, we write the 
value of a on x as a[x], rather than a(x). 



Wc define the adjoint 
: V\ — > V 2 via 



: V 2 * — > V* of a linear map 



(a)[x) = a[(j)(x)} 



(2) 



for all a G V* ,x G V. 

An inner product, written (,), on V distinguishes a 
particular isomorphism £ between V and V*, defined by 
requiring that ((x) G V* satisfy C(x)[y] = (x,y), for all 
y G V. From now on, we will assume such a choice 
of inner product (there will be a natural choice in our 
applications), and canonically identify V with V* via £. 
Thus we will view the dual cone C* as {y G V : (y, x) > 
0,Va; G C}, and the adjoint of (j) : V\ — » V2 as ft : V2 — ► 
V\ defined via 



(B,<j>(A))=(ft(B),A) 



(3) 



for all A e Vi,B e V 2 . 

We say a linear map <j> : V\ — > V2 is Ci-to-C2 positive, 
for cones C\ C Vi, C2 C V2, if 4>{Ci) C C 2 . Under either 
definition of 0^ the following proposition is easily (but 
instructively) verified. 

Proposition 1 If ^{Cx) C C 2 tften ^(C|) C Cf. 



The positive semidefinite cone PSD(d) in the real lin- 
ear space of Hcrmitian dx d matrices, is the set of matri- 
ces M such that x^Mx > for all x G C d . It is self-dual; 
if we use the trace inner product (X,Y) := tr XY to 
identify V* with V, it is not only isomorphic to, but 
equal to, its dual. We will denote by ">r" the partial or- 
der induced by this cone, and often write M y for the 
equivalent M G PSD(d). 

We will have several occasions to use the following 
proposition, which follows from the fact that for normal 
(including Hermitian) matrices, A, ||A||oo is the largest 
modulus of an eigenvalue of A). 

Proposition 2 Let A be Hermitian. Then I + A >z is 
equivalent to HAjjoo < 1. 



III. SEPARABLE CONES 

Let us consider an m-partite unnormalizcd density ma- 
trix (i.e. just positive semidefinite) 

p : H t ® H 2 ® ... ® H m — ► i?i ® H 2 ® ... ® H m 

Let dim(Hi) = di, 1 < i < m. Then any such 

P = {p(h,i2, ■-,i m ;h,j2, -,3m)}, 1 < ik,jk <d k ;l <k < m}. 

Let us block-partition p with respect to the first index : 



P ■= 



/pM p 

P 2 ' 1 P 



1:2 
1,2 



(4) 



^ p^i^ 1 p^i. 2 pdi,di ^ 



where the blocks 
P 

and 



: H 2 



Ho 



H„ 



P l '° = {p(i,h, ...,i m ;j 1 j 2 , ...,j m }, 1 < ik,jk < d k ;2 < k < to}. 



Definition 2 The linear space of N x N complex ma- 
trices is denoted M(N), the linear space over the reals 
of N x N real matrices is denoted Mat(N) , the lin- 
ear space of real symmetric N x N matrices is denoted 
RSym(N), the linear space over reals of N x N complex 
hermitian matrices is denoted as Her(N). The space of 
complex block matrices, K blocks by K blocks, with blocks 
in M(N), is denoted Block(K,N) 
Consider cones d C M(di), 1 < i < to. A matrix 



p : Hi ® H 2 1 



H!®H 2 - 



(i.e. p G M(d\d 2 ...d m ) ) is called (C\ ® C 2 ® ... ® C m )- 
separable if it belongs to the cone generated by the set 
{A x <x> A 2 <g> ... ® A m : A l e d, 1 < i < to} . 
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This is trivially equivalent to the recursive definition: 
Sep(C±, C 2 , C m ) is the cone generated by the pairs 
A 1 <g> B with Ai G d,B G Sep(C 2 , -,C m ). 
Examples: 1. Let for all 1 < i < m the cone Ci 
be the cone of positive semidcfinite matrices, denoted 
by PSD(di). In this case the definition of (Ci <g> C 2 (E> 
• ■■ <8> Cm )-separability is equivalent to the standard no- 
tion of separability of multiparty unnormalized density 
matrices. We will denote the corresponding cone of 
separable multiparty unnormalized density matrices as 
Sep(d 1 ,d 2 , ...,d m ). 

2. Let for all 1 < i < to the cone Ci be the 
cone of positive semidefinite matrices with real entries. 
We call the the corresponding cone the cone of real- 
separable multiparty density matrices and denote it by 



RSep{d\,d 2 , —,d m ). 



: d m ) : 



We now recursively define a subspace RLin(di 
which we show is the minimal linear subspace (over the 
reals) of the symmetric matrices RSym(di, d rn ) that 
contains the real separable cone RSep{d\ 1 d m ). 

Definition 3 p G RLin(di, d 2l d m ) iff in the block- 
partition^ p l ' J G RLin{d 2 , ...,d m ) and p 1,3 = p 3 ' 1 ^- < 
ij j < di); RLin(d) — RSym(d). 

It is easy to prove that p : Hi ® H% ® ... ® H m — > 
Hi ® Hi ® H m is real-separable iff p is separable and 
p G RLin(di, d 2l d m ). In fact, we will show 

Proposition 3 RLin{di,d 2 , ...,d m ) is the minimal lin- 
ear subspace (over the reals) of RSym(did 2 ...d m ) which 
contains RSep{d\ 1 d 2 , d m ). 

Proof: It is clear that RLin(di....d m ) is a subspace of 
("<") RSym(di...d m ). To be explicit, symmetry means 
p(ii,...,i m ,ji,...,j m ) = p(ji,...,j m ,ii,...,i m ); this fol- 
lows from the definition of RLin and RLin(d 2 ...d m Ys be- 
ing a subspace of RSym(d 2 , d m ). That establishes our 
induction step; the base case RLin(d m ) — RSym(d m ) 
clearly also holds. 

Suppose X G RSep(di, ...,d m -i). That is, X = 
J2 k A k (g>B k , A k G RSep(di, ...,d m -i),B k G RSym(d m ). 
By the induction hypothesis, A k G RLin{d\, d m —i). 
Block-partitioning with respect to the second system, 

These blocks are in RLin(di...d m —x) because A k are; 
consequently X G RLin(d\, ...,d m ). The base case is 
trivial: RSep(di) C (RSep(di) n RLin{di)) holds with 
equality because RSep(di) = PSD{d{) and RLin(di) = 
RSym(di) . 

For the opposite direction, let X G 
RLin(di, ...d m ), Sep(d\, rf m ). By separability, 



(6) 



where A k G PSD(d!),B k e PSD(d 2 ),Z k G PSD(d m ) 
(and no restriction to to = 26 is intended!). Let 



A k = A\ + iA\ with A\ real symmetric, A\ real skew- 
symmetric, and similarly for B. Substituting these in © 
and keeping only terms with an even number of imag- 
inary factors (since X G RLin), and block-partitioning 
the matrix according to the first subsystem, each block 
has the form: 



J2Al(i,j)R k +J2Al(i,j)S k 



(7) 



Thus X = Xi + X 2 , where the first term is block- 
symmetric, the second block-skew-symmetric. This sec- 
ond term must therefore be zero. By the recursive def- 
inition of RLin (and Sep) we have that, for each fixed 
value of X lJ must be block-symmetric when parti- 
tioned according to the second ( "B" ) system. This block 
is 



YAl^j)Bl(m,n)C k 



k 



+ yAi(i,j)Bi(rn,n)C« 



k 



■ ■■z k 



...Z" 



(8) 



and again only the first component is nonzero. Proceed- 
ing thus through all the subsystems, all terms with a 
skew-symmetric factor must be zero and we have: 



X = A\ ® B\ 



^k > 



(9) 



with each of A k , B kl . . . Z k real symmetric and positive 
semidefinite, i.e. X G RSep. 

m 

The next lemma gives a simple but very useful crite- 
rion for P5 , Z?(di)(g)C(d2)-separability for any cone C(d 2 ) 
of Hermitian matrices. It is a mild generalization of the 
necessary and sufficient criterion (cf. p|, 0|) for ordi- 
nary (C(d 2 ) = PSD(d 2 )) bipartite separability, that ev- 
ery positive linear map, applied to one subsystem of the 
bipartite system (i.e. to every block of its block density 
matrix) gives a positive semidefinite matrix. 

Definition 4 A linear operator <fi : M(d 2 ) — > M(N) is 
called C(d 2 ) -positive if (f>(C(d 2 )) C PSD(N). If X is a 
block matrix as in O), Xij G M(d 2 ) and <ft ■ M(d 2 ) — ► 
M(N) is a linear operator then we define 



4>(X) 



( HX 1 ' 1 ) 

4>{x^) 



<M* 2 < 2 ) 



4>{x^) \ 

(j>{X 2 ^) 



(10) 



Lemma 1 Suppose that the cone C(d 2 ) C Her(d 2 ) C 
M(d 2 ) . Then X is PSD{d x ) <g> C(d 2 ) -separable iff 
4>{X) y (i.e. is positive semidefinite) for all C(d 2 )- 
positive linear operators <f> : M(d 2 ) — > M(di). 

The proof uses the following proposition, which gener- 
alizes the duality between positive linear maps and sep- 
arable states. 
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Proposition 4 For Hermitian matrices M G 
Block(d%, c# 2 ) , the following are equivalent: 

1. tr MZ > for all P5D(di) ® C(d 2 )- separable 
matrices Z. 

2. M % i = x( e i e j) f or some \ such that 
X (PSD(d 1 )) C C(da)*. 

Proof of Proposition |4] Item 1 of the proposition is 
equivalent to the same statement with Z of the form 
xx> <X> Y, with xx^ a rank-one matrix in PSD(d\) and 
y G C(d 2 ). This is equivalent to: 



V x,x*tr (M y F) =tr N 



^^M'jy^o. (11) 



Since this holds for any Y G C(d2), this says 
J2i XiX*M tj G C(d 2 )*. Define a linear map T M ■ 
M{di) -> M(d 2 ) via T M (e 4 e]) = M y . Then for any 
x G C d \ T(xx^) = Zi j x i x* j T M (e i e]) = V,, r,.r; ,\/'-' , 
which we have just shown is in C(d2)* ■ Since xx* gen- 
erate PSD(di), this is equivalent to saying Tm takes 
PSD(d x ) into C(d 2 )*, so T m is the desired x- The other 
direction (2 =>■ 1) is similar (most of the steps above were 
equivalences). ■ 

We remark that essentially the same argument estab- 
lishes a similar statement for arbitrary pairs of cones of 
Hermitian matrices, in which PSD(d\) is replaced by a 
cone C(d\) in item 1, and by C(d±)* in item 2. 
Proof of Lemma 1: "Only if" is trivial: PSD(d 1 ) ® 
C(G? 2 )-separability of X means X — Yli^-i ® Bi> with 
Ai G PSD(di),Bi G C(d 2 ), hence 0(A) = Y,i A i® 
4>{Bi); since we assumed 0(C(ci 2 )) C PSD(d\), we have 
0(5*) G PSD{dx), so 0(A) G 5ep(di,di). 

For "if," note that 0(X) >: says that for any positive 
semidefinite B G Block(N, K), 



5^tr B ij <f)(X ij ) >0. 



(12) 



By the definition of dual cone (and the self-duality of 
PSD(dx)) it is easily seen that (f>^ (PSD(d 1 j) C C(d 2 )*. 
Now, (|12|l is equivalent to 



^ tr ^(B«)pi = tr 0t(5)X > 



(13) 



Letting B y = e,*ej, so B is the block matrix of a positive 
semidefinite rank-one state (specifically, the unnormal- 
ized maximally entangled state xx< with x = J^i <S) e^), 
we have that the matrix <ft(B) satisfies condition 2 of 
the Proposition; as ranges over all PSD(di)-to-C(d2)*- 
positive maps, (f)^(B) ranges over all such matrices so by 
Proposition E| A is PSD(d 1 ) <£> C(d 2 )-separable. ■ 
The next proposition will allow us to extend the (ex- 
act) bipartite result from 3] to multiparty systems. The 
bipartite result was that everything in the ball B(N, 1) := 
{I + A : ||A|| 2 < 1} of size 1 in Frobenius norm around 



the identity operator is separable; therefore, so is every- 
thing in the cone, which we call G(di<i 2 , 1), generated by 
that ball, for a bipartite system with subsytems of di- 
mensions di 1 d 2 - We define a slight generalization of this 
cone: 

Definition 5 Let G(N, a) C Her(N) C M(N) be the 
cone generated by hermitian N x N matrices of the form 
{/ + A:||A|| a =:(tr(AAt)i <a}. 

A sufficient condition for tripartite separability of X is 
clearly that it belong to the cone generated by Ai® Zi, 
where Ai G PSD(d\) and Zi G G(<i 2 <i3,l); this can be 
used to derive a tripartite sufficient condition for separa- 
bility in terms of Frobenius norm. Letting this tripartite 
2-norm ball generate a cone of separable tripartite states, 
similar reasoning gives a ball of 4-partite states, and so 
on. The key to the induction step is the following propo- 
sition. 

Proposition 5 7/0: M(N) — ► M(K) is a G(N, ex- 
positive linear operator (i.e. 0(A) >z for all X G 
G(N, a) and 0(1) = I G M(K) then 

1. H^JSQIloo < a^ 1 !!^ || 2 for all hermitian X G 
M(N). 

2. ||0(F)|U < a-V2||F|| 2 for all Y G M(N). 

Proof: All G(N, a)-positive satisfy 0(7 + A) > for 
all A such that ||A|| 2 < a, which when (f)(1) — I gives 
I + 0(A) >z 0; by Proposition [21 this is equivalent to 
1 10(A) I loo < 1, establishing item 1. Item 2. uses item 1 
and the following: 

If a linear operator : M(N) — > M(K) satisfies 
(f>(Her(N)) C Her(K) and ||0(Z)||oo < ||A|| 2 for all 
hermitian Z G M(N) then ||0(F)||oo < \/2||Y|| 2 for all 
Y G M(N). 

To show this, let Y = A + iB with A, B Hermitian. 
Now, 0(F) = 0(A) +i(f)(B), so ||0Cn||oc < ||0(^)IU + 
110(5)1100. This is less than or equal to ||A|| 2 + ||73|| 2 
by assumption (A, B being Hermitian). The conclusion 
follows by the square root of the elementary inequality 
(x + y) 2 < 2(x 2 + y 2 ) (obtained from 2xy < x 2 + y 2 , 
which comes from (x — y) 2 > 0). ■ 

The following example shows that the extra v2 factor 
is the best possible in this statement. 

Consider : M (2) — ► M(2), 0(A) = A(l, 1)A 1 + 
A(2, 2)A 2 ; where A\, A2 are real symmetric anticommut- 
ing unitary matrices: 



At 



Ao 



1 

-1 



1 

1 



Notice that for real a, b we have that aAi + bA 2 = 
(a 2 + b 2 )zU for some real symmetric unitary U. Thus 
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(j){Her{2)) C Her(2) and \\(t>{Z)\\oo < \\Z\\ 2 for all her- 
mitian Z £ M(2). Consider the following (nonhermitian) 
matrix : 



Y 



1 

i 



Then = 2 and 11^)11^ = \\HY)\\ 2 2 = 4, since 

Det{(j>{Y)) = 0. ■ 
Problem: Is the extra \[2 factor still the best possible 
if, additionally, (/>(!)_= H 

Proposition 2 of Q is a similar contraction result with 
constant 1 rather than y2 on all matrices, not just Her- 
mitian ones, for the usual positive maps; the proofs used 
there do not work for the different notion of positivity 
used here. 

Now everything is ready for our attack on multipartite 
separability. 



Theorem 1 Let Hi,H 2 have dimensions ni,n 2 . 

unnormalized density matrix p : Hi (g> H 2 ► H\ 

satisfies the inequality \\p — I\ 
PSD(yii) <8> G(n 2 , a)-separable. 



< a/V2 then it 



If an 

® H 2 



Proof: 

Let p = 1 + A A Hermitian; by Lemma 1, we are look- 
ing for a bound on ||A||2 that ensures, for any G(n 2 ,a)- 
to-PSD(ni)-positive linear operator, that <j)(I + A) y 0. 
4>(I) = I, so 4>{I + A) = I + 4>(A); | |0(A) | U < 1 will en- 
sure this (cf. Proposition I2J). The argument establishing 
that H^A^loo < 1 is essentially identical to the proof of 
the main theorem of Q , except that because <ft is not an 
ordinary positive map we must use the weaker contrac- 
tion bound of Proposition [3 with its v2 factor, in place 
of Q's result with a factor 1. 



||0(A)||^< \\A\H < 
where A := [a^], a,j := ||</>(A y )||o< 



(14) 



(15) 



(The first inequality is because the operator norm of a 
block matrix is bounded above by that of the matrix 
whose elements are the norms of the blocks, and the sec- 
ond is because the Frobenius norm is an upper bound 
to the operator norm.) ||0(A y )||^ < 2a- 2 ||A^|| O0 by 
Proposition and it is an elementary norm inequality 
that IIA^Hoo < IIA^Ma. So 



||0(A)||L<2 a - 2 £||0(A^||2 

ij 

<2a- 2 ^||A^||^2a- 2 ||A|| 2 

ij 

Thus if ||A|| 2 < a/V2, \\4>(&)\\oc < 1. 



(16) 



Corollary 1 // an m-partite unnormalized density ma- 
trix p : H± eg) • • • (8) H m — * Hi ® • • • (g) H m satisifes 
\\p- I\\ 2 < l/(2 m / 2 - 1 ) then it is separable. 

Proof: The main result of Q is that G(did 2 , 1) C 
Sep(di,d 2 ). This is the base case for an induc- 
tion on the number of subsystems. For the induc- 
tion step, fix m > 2 and suppose as our induc- 
tion hypothesis the Corollary holds for m — 1, i.e. 
G(di,...,d m _i,2-« OT - 1 >/ 2 - 1 )) C Sep(di, ...d{). Theorem 
□ tells us GH,...^-"™- 1 )/ 2 - 1 )/^ = 2-(™/ 2 - 1 )) 
is P5D(d m )(8G'(di...d m -i,2-( m - 1 )/ 2 - 1 ))-separable, and 
therefore (by the induction hypothesis and the recursive 
definition of separability) separable. ■ 
If we had an analogue to Theorem^ with G{n\, a) in- 
stead of PSD(ni), and some constant a replacing 1/V2, 
then we could get one over a polynomial instead of an 
exponential in Corollary □ by recursively dividing sys- 
tems into subsystems of more or less equal size, since 
this involves a logarithmic number of partitionings com- 
pared to splitting off one system at a time. The first step 
toward such a theorem would be to apply the characteri- 
zation of C(di) ® ^(^-separability analogous to Propo- 
sition □ (and discussed after that proposition above) to 
G(di, a), G(d 2 , a); the fact that neither of these cones is 
self-dual has so far proved an obstacle to our getting use- 
ful results along these lines. 



Theorem 2 Consider p £ RLin(d\ 
I\\ 2 < 1 then 



If\\P 



Pi, a, 



> 



(17) 



where for all i pi is a real positive semidefinite di x d\ 
matrix; pi € RLin(d 2 , d m ) and \\I — Pi\\ 2 < 1. 

Proof: The proof goes essentially like that of Theorem 
□ except that the blocks A 1 - 7 are Hermitian by Proposi- 
tion |3J Consequently we may use item 2. rather than 
item 1. of Proposition [S] and obtain the larger radius 
ball. ■ 

Corollary 2 If p G RLin(d±, ...,d m ) and \\p — I\\ 2 < 
1 then p is real- separable. In other words the maximal 
separable ball in RLin(d±, ...,d m ) around the identity I 
has radius 1 . 

The next Proposition is immediate from results of , 
derived using "scaling," i.e., considering all ways of writ- 
ing a matrix p as a positive scalar times the sum of the 
identity and a Hermitian perturbation, and minimizing 
the 2-norm of the perturbation). 

Proposition 6 Define p(p) as the maximum of ||A||2 
over all A such that there exists an a > for which p = 
a(/+A). Let p be a normalized (tr p = I) density matrix. 
Then the following three statements are equivalent: 
1. p(p) < a. 
2. 
3. 



tr p 2 < l/(d-a 2 ). 

\p-I/d\\ 2 < a/^d(d-a 2 ). 
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Using this Proposition, Theorem ^ has (via Corollary 
^) the following corollary: 

Corollary 3 If an m-partite normalized (i.e. unit trace) 
density matrix p : H\ • • • (8 H m — ► Hi (8 • ■ • (8> -ff m 
satisifes \\p — I / d\\ 2 < 2m/ i_i d , where d = dim{H 1 (g> 
■ ■ ■ ® H m ), then it is separable. 

(The proposition actually gives the (negligibly) tighter 
statement with 2 m / 2 ~ 1 y/d(d - 2"( m ~ 2 )) in the denomi- 
nator.) 

IV. DISCUSSION 

In many interesting experimental or theoretical situa- 
tions, the system is in a "pseudopure state" : a mixture 
of the uniform density matrix with some pure state it: 

p e ^:=€Tr + (l-e)I/d, (18) 

where d — d%, ..d m is the total dimension of the sys- 
tem. For example, consider nuclear magnetic resonance 
(NMR) quantum information-processing (QIP), where 
d = 2 (the Hilbert space of a nuclear spin), and m is the 
number of spins addressed in the molecule being used. As 
discussed in more detail below, the initialization proce- 
dures standard in most NMRQIP implementations pre- 
pare pseudopure states. 
Write 

Pe , n = (l/d)J + £(tt - I/d) . (19) 

Since ||e(7f — I/d)\\2 = e\J ^p, by Corollary 3, this is 
separable if 

e < 2-("*/a-i)/^d(d - i) , (20) 

For m Z?-dimensional systems (so d = D m ), this implies 
the (negligibly loosened) bound 

e < 2~ i - m/2 - x) /D m . (21) 

This is an exponential improvement over the result in 
(the qubit case is in Q) of e < 1/(1 + D 2m_1 ). For 
m qubits, for example, our resultgoes asymptotically as 
2-((3/2)m-l) j versus 2-(2"»-l) m Another comparison 

is with our earlier bound of 

e<l/{D m -l) (22) 

guaranteeing separability for m D-dimensional systems 
with respect to every bipartition Q . It is interesting that 
this is exponentially larger than the present bound guar- 
anteeing multipartite separability, although we do not 
know that a tight multipartite bound would still exhibit 
this exponential separation. 

In liquid-state NMR at high temperature T, the sam- 
ple is placed in a high DC magnetic field. Each spin 



is in a highly mixed thermal state. It is not maximally 
mixed because of the energy splitting between the higher- 
energy state in which the spin is aligned with the mag- 
netic field, and the higher-energy one in which it is anti- 
aligned. This gives probabilities for those states propor- 
tional to the Boltzmann factors e ±l3,lB , where (3 = 1/kT 
with k Boltzmann's constant, fi the magnetic moment 
of the nuclear spin, B the external field strength. For 
realistic high-T liquid NMR values of T = 300 Kelvin, 
B = 11 Tcsla, (3pB » 3.746 x 10" 5 < 1. Calling this r), 
e ±v fa 1 ± 77, so the probabilities are p-f ~ (1 — 77) /2, ~ 
(1 +rj)/2, where ] /[ denote alignment/anti-alignment of 
the spin with the field. With independent, distinguish- 
able nuclear spins, Maxwell-Boltzmann statistics give the 
highest-probability pure state, with all m spins up (field- 
aligned), probability about (1 + r7) m /2 m « (1 + mn)j2 m . 
Standard pseudopure-state preparation creates a mixture 
of this most probable pure state and the maximally mixed 
state, by applying a randomly chosen unitary from the 
group of unitaries fixing the all-spins-aligned state. (U 
could be chosen uniformly (i.e. with Haar measure on 
this group), but efficient randomization procedures may 
draw from carefully chosen finite sets of such unitaries 
6].) Thus, we get a mixture 

(1 - e)I/2 m + e| T ••• TXT (23) 

with 

e = W i/2 m . (24) 

With 7} w 3.746 x 10" 5 . this implies that below about 23 
qubits, NMR pseudopure states are all separable, com- 
pared to the « 13 qubits one gets from the bound in 

Even without entanglement considerations, it is clear 
that pseudopure-state NMR quantum computing will 
not give asymptotic gains over classical computing, be- 
cause of the exponentially increasing signal-to-noise ra- 
tio from As pointed out in Q, it does not fol- 
low that no application of liquid-state NMR can have 
better-than-classical asymptotic performance: NMRQIP 
is not limited to pseudopure-state initialization. It is not 
known that NMRQIP with other initialization schemes, 
such as those that involve preparing a fixed number of 
pseudopure qubits as the total number of spin qubits 
grows, can be efficiently classically simulated, and even 
with one pure qubit interesting things can be efficiently 
done for which no efficient classical algorithm is cur- 
rently known 7]. Another non-pseudopure initializa- 
tion scheme is the Schulman-Vazirani algorithmic cool- 
ing procedure Q , which essentially uses an efficient (and 
NMR-implcmcntable) compression algorithm to convert 
the thermal state with entropy S, for m nuclear spins, 
into log S maximally mixed qubits and m — log S pure 
ones. Their work shows that the theoretical model de- 
rived from NMR with an initial thermal state is as pow- 
erful as standard quantum computation. Though the 
overhead required is polynomial, the space overhead is 
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large enough to be impractical given the small number 
of qubits available in liquid-state NMR. But algorithmic 
cooling is certainly relevant in principle to the asymptotic 
power of an implementation, and could be practically rel- 
evant to a high temperature bulk QIP implementation 
that was expected to be sufficiently scalable that asymp- 
totic considerations are relevant. 

However, there is still the interesting possibility that 
one may produce an entangled overall density matrix 
(and not just a mixture of the maximally mixed state 
with an entangled state) via pseudopure-state NMRQIP. 
The results herein increase the number of qubits known 
to be required before this may be possible, although, 
since we have not shown that the bounds herein are tight, 
with our assumed rj even at 23 qubits there is no guaran- 
tee one can prepare an entangled pseudopure state. By 
contrast, from (|22[) and 124(1 one needs m = qubits 
(about 26, 700 for our rf) to have any hope of obtaining a 
pseudopure state that is not bipartite-separable with re- 
spect to a partition of the qubits into two sets. Again, we 
remind the reader that non-pseudopure protocols could 
conceivably give such a "biseparable" state with far fewer 
qubits; in the Appendix we discuss some implications of 
our results for this possibility. 

In conclusion we have derived an upper bound, ex- 
ponentially better than those already known, on the 
Frobenius-norm radius of a ball of separable matrices 
around the identity matrix. The bound has implications 
for the minimum polarization needed for bulk quantum- 
information-proccssing protocols initialized by preparing 
pseudopure states from a thermal state by via averaging, 
to produce entanglement: known lower bounds on this 
minimal polarization are exponentially increased by our 
results. We stress, however, that this is just one applica- 
tion of a general, computationally simple sufficient crite- 
rion for multipartite entanglement, applicable to states 
of any form. Its geometric nature should make it use- 
ful in many applications, both theoretical and practical. 
The question of whether this bound is tight, or whether 
there is a larger separable Frobenius norm ball around 
the identity, remains open. 



APPENDIX A: ENTANGLEMENT AND 
THERMAL INITIAL STATES IN NMR 

Schulman and Vazirani's algorithmic cooling protocol 
shows that it is, in theory, possible to prepare any entan- 
gled state from sufficiently many thermal NMR qubits. 
The question of just how many qubits are required by 
means possibly simpler than algorithmic cooling is also 
of interest. One can gain some information about this 
using our results, by applying Corollary 3 to the initial 
thermal density matrix of an NMR system. This matrix, 
which is approximately 







2 

±=2 



(with each qubit expressed in the | t), | I) basis), has 

(A2) 



\\p-i/d\\i= [i+ ir 



This should be compared to Corollary 3's separability 
condition, which guarantees separability if this squared 
distance is below 2- ( - 3m - 2 K The comparison gives that 
for m qubits, the thermal state (and any state reachable 
from it by unitary transformation) is separable if 



r) < rrT 1/2 2 r 



(A3) 



For the same experimental conditions considered above, 
14 qubits are required before this bound is exceeded 
(rather than the 23 for the pseudopure state prepared 
from this thermal state). 

If one wants the possibility of bipartite entanglement 
with respect to some partition of the qubits into two sets, 
it is necessary to beat the bound (closely related to Q\22\i ) 



\\p-I/d\\l< 



1 



d{d-l) 



(A4) 



(Al) 



For the thermal state, this gives rj < m _1 / 2 2 - " 1 / 2 . With 
i] = 3.746 x 10~ 5 as before, the bound is not surpassed 
until 25 qubits. Although this comes from an essen- 
tially tight bound, that does not imply that entangle- 
ment can be achieved through computation starting with 
this initial state. Although the thermal state has the 
same magnitude perturbation as some entangled state, 
the latter will have different eigenvalues, so unitary ma- 
nipulation will not get us there, and it is an interest- 
ing question what we can achieve along these lines using 
NMR-implementable nonunitary manipulations (which 
may sometimes require extra thermal ancilla bits that 
should be counted as resources). To understand when 
unitary manipulations might be guaranteed to give us 
entanglement, a promising approach is to look for con- 
ditions on the spectrum of a matrix sufficient for an en- 
tangled state with that spectrum to exist. We conjec- 
ture that the problem of determining, from a spectrum, 
whether or not an entangled state with that spectrum 
exists, is NP-hard in terms of an appropriate measure 
of problem size; this does not rule out easier-to-evaluate 
sufficient conditions, perhaps obtained from relaxations 
of the above problem. Theorem 4 of Q gives some infor- 
mation on spectra sufficient to guarantee entanglement. 

For the pseudopure fraction e there is an upper bound 
of 2/(2 + 2 m ) from Q nearly matching the lower bound 
(12211 . Comparison to the expression e = i]m/2 m for pseu- 
dopure polarization suggests that, at some large num- 
ber of qubits m > 2/?7, even pseudopure protocols will 
exceed this bound. With n ~ 3.746 x 10 -5 this gives 
about 53,400 qubits. This is far beyond the range gen- 
erally viewed as relevant for liquid state NMR, and im- 
proved polarization is unlikely to bring it into this range. 
Since the state demonstrating entanglement at the up- 
per bound is not in general pseudopure, exceeding this 



8 



bound is still no guarantee we can get entanglement. 
Also, this is well into the range where algorithmic cool- 
ing could produce much stronger entanglement. Still, one 
can imagine that in other bulk QIP implementations the 
balance between the difficulty of implementing complex 
unitaries (relatively easy in NMR) and the difficulty of 
preparing large numbers of thermal qubits (apparently 
relatively hard in NMR) could be different, and entan- 
glement generation by simple manipulations on thermal 



states, perhaps even by pseudopure state preparation, 
might be promising. 
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